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Abstract We use all-optical methods to produce a highly-degenerate Fermi gas of 
spin-1 /2 ^Li atoms. A magnetic field tunes the gas near a collisional (Feshbach) 
resonance, producing strong interactions between spin-up and spin-down atoms. 
This atomic gas is a paradigm for strong interactions in nature, and provides tests 
of current quantum many-body calculational methods for diverse systems, includ- 
ing very high temperature superconductors, nuclear matter in neutron stars, and 
the quark-gluon plasma of the Big Bang. We have measured properties of a breath- 
ing mode over a wide range of temperatures. At temperatures both below and well 
above the superfluid transition, the frequency of the mode is nearly constant and 
very close to the hydrodynamic value. However, explaining both the frequency 
and the damping rate in the normal collisional regime has not been achieved. Our 
measurements of the damping rate as a function of the energy of the gas are used 
to estimate an upper bound on the viscosity. Using our new measurements of the 
entropy of the gas, we estimate the ratio of the shear viscosity to the entropy den- 
sity, and compare the result with the lower bound for quantum viscosity recently 
predicted using string theory methods. 
PACS numbers; 03.75.Ss, 32.80.Pj. 



1 Introduction 

Optically-trapped, strongly-interacting atomic Fermi gases''^ provide a unique 
laboratory for testing nonperturbative many-body theories in a variety of fields, 
from neutron stars and nuclear matter -^'^'^ to quark-gluon plasmas^ and high tem- 
perature superconductors^. In contrast to other Fermi systems, in atomic gases, 
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one may tune at will the interaction s energy— li^, and spin population s '^i'"^ . 
The strong interactions are produced using a Fano - Feshbach resonanc e ii'^i^^ , 
where the scattering length for zero-energy s-wave collisions is large compared to 
the interparticle spacing. 

In a strongly-interacting 50-50 mixture of spin-up and spin-down atoms, evi- 
dence of new, high-temperature superfluidity has been found in the studies of col- 
lective dynamic s and confirmed by the observation of quantized vortex 
lattices^". Microscopic properties of the superfluid have been probed by detecting 
pairs of Fermi atoms in projection experiments radiofrequency spectroscopy 
experiments^^'^'*, and optical spectroscopy measurements of the order parame- 
ter—. Measurements of the heat capacity collective damping^^, and entropy^- 
near resonance reveal transitions in behavior close to the temperature predicted 
for the onset of Fermi superfluidity. Despite the microkelvin range of tempera- 
tures, the phenomenon is referred to as "high-temperature superfluidity" because 
the transition temperature is a large fraction (ci 30%) of the Fermi temperature. 

Measurements of the physical properties of strongly interacting Fermi atoms 
are not limited to atomic physics and can be extended to other Fermi systems: 
Fermi gases near a broad Feshbach resonance exhibit universal interaction s 
and universal thermodynamics i. e. the properties of the gas do not depend 
on the details of the interaction potentials (as long as the potential range is small) 
and are identical to those of other resonantly-interacting Fermi systems. For exam- 
ple, in a uniform strongly interacting gas, the ground-state energy is a universal 
fraction, denoted 1 + j3, of the energy of a noninteracting gas at the same den- 
sity The universal energy relationship was originally explored theoretically in 

the context of nuclear matter^^ and has later been measured using ultracold Fermi 
atom^*2i2J±B3i. 

The strongly-interacting atomic gas behaves as a fluid over a wide range of 
temperatures In the lower temperature range, this hydrodynamic behavior is 
explained by superfluidity At temperatures well above the superfluid transi- 
tion, the observed hydrodynamic frequency and the measured damping rate are 
not consistent with a model of a collisional normal gas 12,33,34 -pj^g microscopic 
mechanism that provides hydrodynamic properties at high temperatures is an open 
question. 

In this paper, we review the hydrodynamic properties of a strongly-interacting 
gas consisting of an equal mixture of spin-up and spin-down atoms. We also 
present model-independent data on hydrodynamic oscillations of the gas, i.e., the 
frequency and damping rate of the radial breathing mode as a function of the en- 
ergy. These results should be useful for future comparison with theoretical models. 

The search for a fluid with the viscosity at the quantum minimum is an- 
other feature of our study of strongly-interacting atomic Fermi gases. On general 
grounds, Kovtun et al.^^, predicted that for any fluid the ratio of shear viscosity 
77 to entropy density i is always > h/AnkB {kg is Boltzmann's constant). For a 
resonantly-interacting atomic gas, we estimate the ratio r\/s from previous mea- 
surements of damping'^ and entropy We discuss whether the atomic gas can 
be used for testing the fundamental r\/s limit. 
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2 Experimental system 

The starting point of our experiments is a highly degenerate, strongly-interacting 
Fermi gas of ^Li. The cloud is prepared using evaporation of an optically-trapped, 
50-50 mixture of spin-up/down states at 840 G, just above the center of a broad 
Feshbach resonance'- To reduce the temperature, we do not employ a 
magnetic sweep from a molecular BEC, in contrast to several other group s ^^i^^i^^ ^ 
Instead, we evaporate directly in the strongly attractive regime, nearly on the Fes- 
hbach resonance ' . 

During the forced evaporation, the depth of the CO2 laser optical trap is re- 
duced to a small fraction (0.002 - 0.0005) of its maximum value (c^i 700 ;UK), 
and then recompressed to the final value, usually 4.6% of the maximum trap 
depth for the breathing mode experiments From the trap frequencies measured 
at 4.6% of full trap depth, and corrected for anharmonicity, we obtain the trap 
aspect ratio X = (0z/C0i_ = 0.045 (O)^ = ^(O^COy) and the mean oscillation fre- 
quency ft) = {(OxCOyO^y^^ =27ix 589(5) Hz. The shape of the trap slightly departs 
from cylindrical symmetry: C0x/(0y = 1.107(0.004). Typically, the total number 
of atoms after cooling is A' = 2.0(0.2) x 10^. The corresponding Fermi energy 
Ef = ksTp and Fermi temperature Tf at the trap center for a noninteracting gas 
are given by 7> = {3NYI^h(b/kB ^ 2.4 jUK at 4.6% of the maximum trap depth. 
For these conditions, the coupling parameter of the strongly-interacting gas at 
B — 840 G is kpa ^ —30.0, where hkp = yJlrnksTp is the Fermi momentum, and 
a — a{B) is the zero-energy scattering length estimated from the measurements of 
Bartenstein et al.— . 

This completes preparation of the gas at nearly the ground state. The energy of 
the gas is increased from the ground state value by abruptly releasing the cloud and 
then recapturing it after a short expansion time them- During the expansion time, 
the total kinetic and interaction energy is conserved. When the trapping potential 
U (x) is reinstated, the potential energy of the expanded gas is larger than that of 
the initially trapped gas, increasing the total energy by a known factor—. After 
waiting for the cloud to reach equilibrium, the sample is ready for subsequent 
measurements. The energy is then directly measured in model-independent way 
from the mean square axial cloud size as described below. 

3 The thermodynamic parameter: Energy or Temperature 

Equilibrium thermodynamic properties of the trapped gas, as well as dynamical 
properties, can be measured as functions of the temperature T or of the total en- 
ergy Efo, . Knowledge of one such variable completely determines the properties 
of a universal system with a known particle number N (provided the local density 
approximation is valid). 

To date, two model-independent thermometry methods have been reported for 
a strongly-interacting gas. The thermometry method of the MIT group is only 
applicable to imbalanced mixtures of spin-up and spin-down atoms and is based 
on fitting the non-interacting edge of the majority cloud by a Thomas-Fermi dis- 
tribution with unconstrained temperature. The method of the Duke group is based 
on measuring the entropy as a function of energy, Stot{Etot)^^ ■ This method has 
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been reported for a balanced mixture of spin-up and spin-down fermions. Noting 
that l/T = dStot / dEtot, it has been proposed that for known energy the tempera- 
ture can be obtained by differentiating a fit to the Stot{Etot) data. In this method, 
however, it is necessary to parametrize the data. 

The other two known thermometry methods are model dependent and rely on 
magnetic field sweeps between the molecular BEC or noninteracting gas regime 
and the strongly interacting regim e . The temperature of the strongly interact- 
ing gas is then estimated from that measured in the BEC or noninteracting gas 
using a theoretical model of the entropy"*". The other thermometry method of the 
Duke groupii. is based on comparing the measured density distribution with an 
approximate model for the density profiles^'. 

A model-independent parametrization of the data is provided by measuring 
the energy of the trapped gas. In the universal regime, the energy per particle E = 
Etot /N can be easily measured because the trapped gas obeys the virial theorem—: 

(f/(x)) = f, (1) 

where U (x) is the harmonic trapping potential and {U (x)) represents average po- 
tential energy per particle. This result is remarkable: Despite the complicated 
many-body strongly interacting ground state, the gas obeys the same virial the- 
orem as a non-interacting Fermi gas. This theorem provides a simple model- 
independent energy measurement, just from the size of the cloud. In the local 
density approximation, with an isotropic pressure, it is easy to show from Eq. [T] 
that the energy per particle in a 50-50 mixture of two spin states in a harmonic 
trap is 

E = 3mO)^{z^), (2) 

where z is the axial direction of the cigar-shaped cloud. The energy can also be 
measured from the radial dimension of the expanded cloud in the hydrodynamic 
regime, where the expansion factor is known ' . 

For the radial breathing mode near the Feshbach resonance, we present both 
the frequency and damping rate data as a function of energy per particle to provide 
model-independent results. 



4 Hydrodynamic flow of a strongly interacting Fermi gas 

Hydrodynamic flow was observed in the very first and simplest experiments with 
a strongly-interacting Fermi gas, where the gas was released from a cigar-shape 
trap ^ . The signature of hydrodynamic flow was the anisotropic expansion of the 
cigar cloud into a disc, i.e. the flow was in the direction of the largest pressure 
gradient, inverting the aspect ratio after expansion. 

In-trap hydrodynamics can be studied by observing the compression (breath- 
ing) mode of the gas ^■'2'^^ '^''^. There are two distinct frequency regimes, for 
a hydrodynamic fluid and for a collisionless normal gas. In a cigar-shape trap 
with our parameters, for the radial mode, these frequencies are ft) = 1.84ft)^ and 
ft) = 2.10ft)^, respectively. In the universal gas, the hydrodynamic frequency has 
been predicted to stay the same at any temperature as long as the flow remains 
isentropic^^. 
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In the experiment, the radial breathing mode is excited by releasing the cloud 
and recapturing the atoms after 25 ^s (for 4.6% of the maximum trap depth). 
After the excitation, we let the cloud oscillate for a variable time tf,oM, at the 
end of which the gas is released and imaged after c^i 1 ms of expansion'^. The 
snap-shots of the oscillating gas are shown in Fig. [T] The measured frequency 




• till 



Fig. 1 Snap-shots of gas after oscillation for a variable time ti,oid, followed by release and ex- 
pansion for 1 ms. Each measurement is destructive. ti,oiii is increasing from left to right. 



(Fig. |2| remains at the hydrodynamic value over the nominal temperature range 



T = 0.12— 1.1 Tf , which corresponds to an energy range from nearly the ground 
state value ~ 0.5 Ef to 3.0Ef. The frequency stays far below the frequency of a 



3 



2.0- 



O 



- ■ - ■ » 



0.5 



1.0 



1.5 



E/E^ 



2.0 



2.5 



3.0 



Fig. 2 Breathing mode in a strongly-interacting gas. Normalized frequency versus the normal- 
ized energy per particle. The superfluid phase transition is located at Ec — Ep. The upper dotted 
line shows the frequency for a noninteracting gas. 



non-interacting gas. No signatures of the superfluid transition are seen in the fre- 
quency dependence, but we can obtain several estimates of the critical energy from 
measurements of other thermodynamic quantities, such as the entropy and heat ca- 
pacity. In these two quantities, a change in behavior is observed at Ec = 0.94Ef 
(T = 0.29 Tf) and E^ = 0.85 Ef (T = 0.21 Tf), respectively. In Figure O we dis- 
play the normalized damping rate 1 /(0±_ T versus the energy. In contrast to the 
frequency, the damping rate exhibits interesting features. AtEc = 1.01 Ef one ob- 
serves a clear change in the damping versus energy dependance: The monotonic 
rise switches to flat dependance, which might be a signature of a phase transition. 

In examining these features, it is appropriate to ask: What are the microscopic 
mechanisms that cause the gas to have hydrodynamic properties? Hydrodynamic 
behavior can appear via at least two mechanisms: (i) superfluidity and (ii) normal 
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Fig. 3 Breathing mode in a strongly-interacting gas. Normalized damping versus the normalized 
energy per particle. 
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dynamics of atomic gas with large number of collisions. Below Ec = 1.01 Ef, the 
reduction in the damping rate as T ^ is consistent with superfluidity, and incon- 
sistent with the scenario of a collisional normal gas '2.17 ^|,ove Ec, however, colli- 
sional dynamics of a normal atomic gas does not completely explain the observed 
behavior. For such a normal system, Bruun and Smith^^ found that hydrodynamic 
behavior is only possible at T < 1 .4 Tp, i.e. only for temperatures where the colli- 
sion rate is large enough. The highest observed damping ratel / t = 0. 12 O)^ occurs 
for E > 2.2 Ef and is consistent with predictions for a transition from hydrody- 
namic to collisionless behavior^^. However, the high damping rate is predicted to 
occur with a higher frequency, closer to 2. 10 (D^ 

It is possible that above the phase transition, the dynamics is significantly af- 
fected by the presence of non-condensed pairs. The step in damping rate at about 
E 2Ef might be due to pair-breaking since at this energy, the trap-averaged 
binding energy becomes smaller than hxo '2. A pseudogap formalism'''^' as well 
as the recent observation of the MIT group suggests that atoms can be paired at 
temperatures well above the phase transition. 



5 Quantum viscosity 

In a strongly interacting system Fermi gas, where the interparticle separation I <x 
1/kf sets the length scale, there is a natural unit of shear viscosity T] which has 
dimension of momentum divided by cross section. The relevant momentum is the 
Fermi momentum, hkf = h/l. The relevant area is determined by the unitarity- 
limited collision cross section 4n/kp l^. Hence, 77 h/P =hn, where n is the 
local total density, so that 

T]=ahn. (3) 

Here, a is generally a dimensionless function of the local reduced temperature 
T/Tp{n), where Tf(n) = tr' (3 n)^/^ /2mkB is the local Fermi temperature. Eq.[3] 
has been discussed by Shuryak^-^. We see that viscosity has a natural quantum 
scale, hn. If the coefficient a is of order unity or smaller, we say that the system 
is in the quantum viscosity regime. 
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Using string theory methods, Kovtun et al., have shown that for a wide class of 
strongly interacting quantum fields, the ratio of the shear viscosity to the entropy 
density has a universal minimum value The entropy density 5 has units of nkg. 
Hence, in the ratio ri/s the density cancels so that ri/s has natural units of h/ks. 
The string theory prediction is 

i 471 kg 

An important question is then how close a strongly interacting Fermi gas 
comes to the minimum quantum viscosity limit. Answering this question requires 
the determination of two physical quantities, shear viscosity and entropy density. 
To estimate the ratio 77 /s, we separately integrate the numerator and denominator 
over the trap volume. Then using / d?y.n = N, where A'^ is the total number of 
atoms and j cPxs = Stot is the total entropy, one obtains, 

s Jd^xs kgS/kB' 

Here, (a) = (1 /N) J d^xna{x) is the trap average of the dimensionless universal 
function a that determines the shear viscosity according to Eq.[3]and S = Stot/N 
is the entropy per particle which has been measured as a function of energy 

We estimate {a){E) from the damping rate 1/t of the radial breathing mode, 
Fig. [3] If the damping rate arises from shear viscosity, then the value of (a) is read- 
ily determined. Using the shear pressure tensor in the hydrodynamic equations 
of motion for the radial breathing mode^^, we easily obtain 

' '-^lay (6) 
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Fig. 4 Quantum viscosity in strongly-interacting Fermi gas. The local shear viscosity takes the 
form T] = ahn. In the figure, (a) is a trap-averaged value of the dimensionless parameter a. 



We determine the energy per particle in a model independent way by exploit- 
ing the virial theorem, as discussed above Using our measured damping ra- 
tios 1 / {t(0±) as a function of energy, we determine (a) {E), as shown in Fig. |4] 



8 



Using the values of («)(£) from Fig.|4]and the measured entropy S{E) from 
Ref.— in Eq.[5] we estimate the ratio of the viscosity to entropy density, as shown 
in Fig. |5] A similar plot has been given previously by Thomas Schafer, based on 
our damping and entropy measurements^^. In the plot, note that Eq :^ 0.5 Bp is 
the ground state energy and Ep is the Fermi energy of a noninteracting gas at 
the trap center. The superfluid transition occurs near E = Ep, above which the 
gas is normal. For comparison, the estimated value for ^He and '^He near the A- 
point i& r\/s 0.7. For a quark-gluon plasma, a current theoretical estimate^ is 
77/^ = 0.16-0.24. 




Fig. 5 Ratio of the shear viscosity r\ to the entropy density s for a strongly interacting Fermi gas 
as a function of energy E, red solid circles. The lower green dotted line shows the string theory 
prediction l/(4;r). The light blue bar shows the estimate for a quark-gluon plasma (QGP)i^, 
while the blue solid bar shows the estimate for ^He and *He, near the X- point. 



While our initial estimates for the r\/s ratio place the strongly interacting 
Fermi gas in the quantum viscosity regime, a caveat for our result is that we have 
not proven that the damping arises from viscosity'^. Other sources of damping 
may contribute, for example Landau damping'*^. Thus, the actual viscosity may 
be lower. Our data may not be determining the minimum viscosity, which may 
make a contribution that is smaller than the total measured damping rate. 
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6 Conclusion 

We have explored the fluid properties of a strongly interacting Fermi gas. Mea- 
surements of the frequency and damping of the radial breathing mode as a func- 
tion of energy provide model-independent data that indicate very low viscosity, 
fluid-like behavior over a wide range of energies covering both the superfluid 
and normal phases. The measured frequency agrees with the macroscopic model 
of tiniversal isentropic hydrodynamics. The unitary collision dynamics that pro- 
vides this fluid-like behavior well above the phase transition is not completely 
understood. Our estimate of the viscosity over entropy density ratio (TJ /s) places 
strongly interacting Fermi gases in the quantum viscosity regime. 
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